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Abstract 


The Global Reference Atmospheric Model as developed anii .e- 
ported in NASA TMX-64871 and 64872 is used along with the revised 
perturbation statistics (January 1975 report on Contract NAS8-30657) 
to evaluate and computer graph various atmospheric statistics along 
a Space Shuttle Reference Mission and Abort Trajectory furnished by 
NASA. The trajectory plots are height vs. ground range, with height 
from ground level to T55 km and ground range along the re-entry tra- 
jectory. 

Cross-sectional plots, height vs. latitude or longitude, are 
also generated for 80® longitude, with heights from 30 km to 90 km 
and latitude from -90° to +90°, and for 45° latitude, with heights 
from 30 km to 90 km and longitudes from 180° E to 180° W. 

The variables plotted are monthly average pressure, density, 
temperature, wind components, and wind speed and standard deviations 

and 99^^ inter-percentile range for each fjf these variables. 
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1 . INTRODUCTION 

Under a previous contract (NAS8-29753) an extended four Dimensional 
(height, latitude, longitude, and month-by-month variation) Global Refer- 
ence Atmospheric Model for surface to orbital altitudes has been developed 
and tested at Georgia Tech. The Global Reference Atmospheric Model com- 
bines: '{) the NASA 4-D worldwide atmospheric model (0 - 25 km) developed 

by Allied Research Associates (NASA CR-2082), 2) the Jacchia 1970 model 

(Smithsonian Astrophysical Observatory Space Report 313) for the thermo- 
sphere above 90 km, and 3) the Groves model (AD 737 794) for latitude 
and longitude dependent stationary perturbations developed at Georgia 
Tech from readings of 10 mb, 2 mb, and 0.4 mb map data and extrapolations 
to 90 km by the method developed for NASA by Northrop (NASA CR-2223). 

In addition to the monthly means, quasi-biennial and random pertur- 
bations in pressure, density, temperature, and winds are evaluated. The 
quasi-biennial amplitudes and phase were determined at Georgia Tech from 
analysis of Meteorological Rocket Network data from 1964 through 1969 
(Justus and Woodrum, NASA CR-2203) and from literature reports of other 
investigators. 

The Global Reference Atmospheric Model as developed and reported in 
NASA TMX-64871 and 64872 is used, along with the revised perturbation 
statistics (January 1975 report on Contract NAS8-30657) to evaluate and 
computer graph various atmospheric statistics along a Space Shuttle Refer- 
ence Mission and Abort Model Trajectory furnished by NASA; The trajectory 
plots are height vs. ground range, with height from ground level to 155 km, 
and ground range along the re-entry trajectory. The 30 km to 90 km height 
range is plotted for January, April, July and October and the 0 km to 30 km 
height range for January and July. The 90 km to 155 km height range is 


plotted for January and July for highi average, and low li^vels of solar 
activity. ^ 

Cross-sectional plots, height vs. latitude or longitude, are also 
developed for 80° longitude, with heights from 30 km to 90 km and latitude 
from -90° to +90°, and for 45° latitude, with heights from 30 km to 90 km 
and longitudes from 180° E to 180° W. These cross-sections were plotted 
for January, April, July, ancl October. The variables plotted are monthly 
average pressure, deni>ity, temperature, wind components, and wind speed 
and standard deviation and 99th inter-percentile range for each of these 
vari ables. 


2 


2. DESCRIPTION OF THE PLOT PROCEDURE 


The plot program which was developed under the present contract will 
generate a plot of any of twelve different variables against a selection 
of X and y axis parameters. It requires as input a data file, which is 
generated by GRAM at appropriate locations; one input card for each variable 
to be plotted; and, if a range plot is to be generated, a range input file, 
which contains the range and corresponding time and is used to label the 
time axis and mark the line of the trajectory. 

The input card for each plot contains the following parameters: 


Parameter 


Definition 


IX 

lY 

IP 

CONT(l) 

C0NT(2) 

LOG 

KON 

IDELTA(7) 

LIMITS(7) 

LEN(7) 

INITIAL(7) 

i 

ISOL 


the parameter for the x axis 

the parameter for the y axis 

the parameters to be plotted 

Contour value for the first plot parameter 

Contour value for the second plot parameter 

the value, in degrees, of the latitude/longitude for a 
height vs. longitude/latitude plot 

specifies whether latitude or longitude is held con- 
stant for a height vs. latitude/longitude plot. LOC 
and KON are used to label the plot 

incremental change along range axis, km 

number of data values along the range axis minus 1 

length of range axis, inches 

initial value of range - used in labeling range axis. 
IDELTA(7), LIMITS(7), LEN(7), and INITIAL(7y are input 
only for ground range plots 

solar activity parameter used to label plot. 1. high 
solar activity, 2. average solar activity, 3. low solar 
acti vi ty 


3 


The following values are permis sable for IX and lY on the input card 
(Note the axis associated with that value): '' 


IX or lY Definition Axi s 


1 

Latitude -90 to +90 

X - axis 

2 

Latitude -90 to +90 

y - axis 

3 

Height 0 to 30 

y - axis 

4 

Height 30 to 90 

y - axis 

5 

Height 90 to 155 

y - axis 

6 

Height 155 to 200 

y - axis 

7 

Ground range "input" to "input" 

X - axis 

8 

Longitude 180 to 0 

X - axis 

9 

Longitude 0 to -180 

X - axis 


The plot parameter, IP, values can be selected from the following list 
Value Definition 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Pressure {% deviation from U.S. Standard) and standard 
deviation of pressure {% relative to GRAM monthly mean) 

Density {% deviation from U.S, Standard) and standard 
deviation of density {% relative to GRAM monthly mean) 

Temperature and standard deviation of temperature (° K) 

Geos trophic monthly mean and standard deviation of east- 
ward wind (m/s) 

Geostrophic monthly mean and standard deviation of 
northward wind (m/s) 

Wind speed and standard deviation of wind speed (m/s) 

•hh 

Upper and lower 99^^ percentile of pressure (% deviation 
from U.S. Standard) 

■m 

Upper and lower 99^^ percentile of density (% deviation 
from U.S. Standard) 

Upper and lower 99^^ percentile of temperature (° K) 

‘hh 

Upper and lower 99 " percentile of eastward wind (m/s) 

XU 

Upper and lower 99 ^ percentile of northward wind (m/s) 

XU 

Upper and lower 99^^ percentile of wind speed (m/s) 







The first parameter for each. IP value appears as a solid line on the 
plot and the second parameter appears as a dashed line. 

The data values used to plot pressure, density, temperature, standard 
deviation of pressure and standard deviation of density are taken directly 
from the GRAM output without change. The standard deviation of temperature 
is changed from % relative to GRAM monthly mean to ° K by: 

Oj = (T)(%y)/100 where csj - standard deviation of tem- 
perature ° K 

r - GRAM monthly mean 

%j - 7a relative to GRAM monthly 
mean, as output by GRAM. 

The wind components and their standard deviation are in m/s as output by 
GRAM. 

The wind speed and standard deviation are calculated as follows: 


( 1 ) 


( 2 ) 


S - / u^ + v^ 
w 


where 


Is = Cl 2 + a/ 

w 


where 


S, , - wind speed 
w 

u - mean eastward wind 
V - mean northward wind 
Tp - standard deviation of wind 


'w 


speed 


0 - standard deviation of east- 

” ward wind 


standard deviation of north- 
ward wind 


Note that the wind speed defined by ( 1 ) is the magnitude of the resultant 
wind and is neither the average, nor the median (50°/) wind speed. 

The 99^^ percentile values for pressure and density are in percent 
deviation from the standard atmosphere as follows: 


/ / 
1 / 
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Zgg - 99^^ percentile value 

(3) Zgg = Z + 2.326 Z - GRAM monthly me^n 

a - standard deviation of GRAM 
monthly mean 

Note that, due to the distinct non-Gaussian nature of the wind speed distri- 
bution* equation (3) does not actually give the 99^" inter-percentile wind 
speeds (1 percentile speeds computed as less than zero by (3) are taken to 
be zero). From Appendix A, equation (3) should give approximately the 99.6^ 
inter-percentile wind speed. 

The percent deviations are obtained by: 


^99 " ^s 


’99 


X 100 whehe Z^ - standard atmosphere value 

%QQ - the respective 99^^ per- 
centile value in percent 
deviation from U.S. Stan- 
dard. 


,th 


The 99 percentile values of temperature are in ° K and tnoSe of the 
wind components are in m/s. They are calculated from equation (3) with mean 
plus standard deviation having the appropriate units. The wind speed per- 
centiles are calculated from (3) with mean and standard deviation as in (1) 
and (2) respectively. 

Cross-sectional plots at 80° longitude are shown in Figures 1 - 12 
for January and in Figures 13 - 24 for April. For July and October plots, 
reverse the signs on the latitude values of the January and April plots 
respectively. All plots are for 30 km to 90 km altitude and -90° to +90° 
latitude. 

Gross-sectional plots at 45° latitude are shown for January (Figures 
25 - 48), April (Figures 49 - 72), July (Figures 73 - 97) and October (Fig- 
ures 97 - 120). The first 12 plots in each month group are for 0° to 180° W 
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longitude and the second 12 are for 180° to 360° W longitude. All plots are 
for 30 km to 90 km altitude. 

The OFT re-entry trajectory provided by NASA is shown in Table 1. 
Trajectory plots are shown for 3 height ranges, 0 km to 30 k?n, 30 km to 

■A 

90 km, and 90 km to 155 km. The 0 km to 30 km height range have plots for 
January (Figures 121 - 132) and July (Figures 175 - 186) while the 30 km to 
90 km range has plots for January (Figures 133 - 144), April (Figures 163 - 
174), July (Figures 187 - 198) and October (Figures 217 - 228). The wind 
plots were not generated for the 90 km to 155 km height range, but different 
levels of solar activity were used to generate January (Figures T45 - 162) 
and July (Figures 199 - 216) plots for this height range. Values used for 
the three levels of solar activity are as follows: (1) high solar 
activity - F-jq -j = 230, F-jq ^ = 230, Ap = 400 , local time = 1400 at 155 km 
point, (2) average solar activity - F-jq ^ = 150, F-jq ^ = 150, Ap = 20.3, 
local time = 0900 at 155 km point, (3) low solar activity F-jq -j - 68, 

F-jq ^ = 68, Ap = 0, local time = 0400 at 155 km point. 

DISCUSSION OF RESULTS 

The figures presented in this report can readily be used to determine 
which seasons would be best or worst for various types of re-entry orbits. 

For example. Figures- 435, 165, 189, arid 219 for the OFT orbit density in Jan- 
uary, April, July and October show the largest density gradients in the 70 - 

I 

90 km altitude range to be in January (from +11 to -9 % U.S. Standard in about 
50 seconds), while the best month is July (from -1 to + 7 5^ U.S. Standard in 
50 seconds). Figures 3 and 15 show that for a polar orbit July would produce 
the strongest density gradients, while April and October would be best, and 
January would be only slightly better than July. 
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Table 1 


OFT re-entry trajectory. Time is in seconds, height 
in kilometers, and latitude and West longitude in 
degrees. The first position shown was used as the 
zero - range position. / 
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APPENDIX A 

As described by 0. E. Smith in NASA TNIX-73319, the probability dis- 

tirbution of wind speeds is a generalized Rayleigh distribution (coming 

from a bi-variate Gaussian distribution of wind components). If we assume 

for convenience that ar>d thpt the correlation <uv> is zero, then 

the X value corresponding to a given probability p is, from equation (12) 

« 

of NASA TMX-73319 

X^ = -2 -en (1 - p) (A-1) 

2 - - 
where x is given in terms of u, v, and a by 

= [(u - u)^ + (V - v)^]/o^ (A-2) 

i.e. X is the non-dimensional radius of the circle enclosing p fraction of 
the wind speed vectors. Since the wind vector (u, v) is the resultant wind, 
with speed V^, the upper speed corresponding to probability p would be 
+ Xa and the lower speed would be - Xa. For p = 99 % x would be 
3.035 or in terms of , where = "^v^ = /T a the upper and 

lower limits on speed would be 

- V + = + 2.146 (A-3) 

I Vi# 

Since the upper and lower percentile wind speeds were computed from equat- 
ion (3) computed by 

V + = + 2.326 0^ (A-4) 

this corresponds to a x value of 3.289 or a probability p of about 99.65^. 

Note that when is negative, zero is used as the lower 99 percentile speed. 
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